General Relativistic effects on the conversion of nuclear to two-flavour quark matter 

in compact stars 
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We investigate the General Relativistic (GR) effects on the conversion from nuclear to two-flavour 
quark matter in compact stars, both static as well as rotating. We find that GR effects lead to 
qualitative differences in rotating stars, indicating the inadequacy of non-relativistic (NR) or even 
Special Relativistic (SR) treatments for these cases. 
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Strange Quark Matter (SQM), consisting of approxi- 
mately equal numbers of up (m), down (d) and strange 
(s) quarks, is the putative true ground state of strong 
interaction [l| , a conjecture supported by model calcula- 
tions for certain ranges of values for strange quark mass 
and strong coupling constant 2]. There have been con- 
certed efforts at confirming the existence of Quark-Gluon 
Plasma (QGP) and SQM, though transiently, in ultra rel- 
ativistic collisions. On the other hand, QGP and SQM 
could naturally occur in the cores of compact stars, where 
central densities of about an order of magnitude higher 
than the nuclear matter saturation density are expected. 
Given the very low strangeness fraction in normal nuclear 
matter (NM), it is almost inevitable that a transition 
from nuclear (hadronic) to quark matter should proceed 
through a conversion to an initial stage of (metastable) 
two flavour quark matter, which should decay to the sta- 
ble SQM. Thus, neutron stars with sufficiently high cen- 
tral densities ought to get converted to strange, or at 
least hybrid, stars. These transitions could have observ- 
able signatures in the form of a jump in the breaking 
index and gamma ray bursts [3, 0] ■ On the other hand a 
full quark star may explain the phenomena of observed 
quasi periodic oscillations [5] . 

There are several plausible scenarios where neutron 
stars could convert to quark stars, through a "seed" of 
external SQM or triggered by the rise in the cen- 
tral density due to a sudden spin-down in older neutron 
stars Several authors have studied the conversion 
of nuclear matter to stran ge m atter under different as- 
sumptions @, i, M M MM M M, M M ElLThey 

have been summarized in a recent work of ours [l9| and 
for the constraint of space, we do not re pea t them here, 
except to mention that Tokareva et al [15| have lately 



modelled the hadron to SQM conversion as a single step 
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process, arguing that the mode of conversion would vary 
with the temperature as well as the value of bag constant. 
Berezhiani et al [l6| . Bombaci et al (l7j and Drago et al 
suggested that the formation of SQM may be de- 
layed, if the deconfinement process takes place through 
a first order transition [20| so that the purely hadronic 
star can spend some time as a metastable object. 

In our recent work [Toj . we have argued that the con- 
version process is a two step process. The first process in- 
volves the deconfinement of nuclear to two-flavour quark 
matter and the conversion process takes some millisec- 
onds to occur. The second process deals with the con- 
version of excess down quarks to strange quarks, which 
occurs via weak interaction, forming a stable SQM and 
the time taken for this process to occur is of the order 
of 100 seconds. This work was restricted to the case of 
a static neutron star and only a SR treatment. To the 
best of our knowledge, GR effects in such processes have 
not yet been addressed in the literature. In this letter, 
we consider GR effects in both static and rotating stars 
and find qualitatively new results. 

As in [101 , we use the Nonlinear Walecka model for the 
nuclear matter euqation of state (EOS). For the sake of 
brevity, we do not repeat the details here. Suffice it to 
say, the star is assumed to be composed of only nucleons. 
The metric describing the structure of the star, is given 
by [21] 



e^-Pr^sir?e{d(l)-ujdtj^ (1) 



-^d9^) 



The four gravitational potentials a,j,p and uj are func- 
tions of 6 and r only. Once these potentials are known, 
one can calculate the observed properties of the star. We 
solve the Einstein's equations for the three potentials 7, p 
and uj, using Green's function technique [1,[2^ and deter- 
mine the fourth potential a from the other three poten- 
tials. The solution of the star is obtained from the 'rns' 
code [23*1, which requires the EOS and a central density 
as inputs and returns various gravitational potentials and 
hydrodynamic parameters as outputs. We solve this code 
for a static as well as for rotating stars with different ve- 
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locities, upto the mass-shedding hmit or the Keplerian 
velocity. The results are, of course, as expected; rotation 
induces a change in the star to an oblate-spheroid shape. 

To study the GR effects on the conversion of nuclear 
to two flavour quark matter, we heuristically assume the 
existence of a combustive phase transition front of in- 
finitesimal thickness, and study the outward propaga- 
tion of the front, using GR hydrodynamical equations of 
motion. First we consider a non-rotating neutron star, 
where the geometry of the problem is one-dimensional 
and then extend our calculation to a rotating neutron 
star, where an extra dimension is needed. Let us assume 
that the conversion front is generated at the center of 
the star, and it propagates outwards through the star 
with a certain velocity, dictated by the initial conditions 
and hydrodynamical equations, converting the nuclear 
matter to two-flavour quark matter. Employing the con- 
servation conditions [2^ and further employing the 
entropy condition 26] , we determine the fLow velocity of 
matter in the two phases, given by 



and 



(P2 -Pl)(£2 

(£2 - ei)(ei +P2) ' 



{P2 -Pi)(ei +P2) 
(62 - ei)(e2 + Pi) ' 



(2) 



(3) 



It is possible to classify the various conversion mecha- 
nisms by comparing the velocities of the respective phases 
with the corresponding velocities of sound, denoted by 
Csi, in these phases. For the conversion to be physically 
possible, velocities should satisfy an additional condition, 
namely, < w| < 1; here we have used natural units 
h, = c = K — 1. If we plot different velocities against 
baryon number, we get curves from which we could cal- 
culate the initial velocity of the front at the centre of 
the star [T9| . Starting with this initial velocity, we in- 
vestigate the evolution of the front with time under GR 
effect; to this end, we first have to derive the appropriate 
continuity and Euler's equations for the metric [eq. (1)]. 
Treating both nuclear and strange matters as ideal fluids, 
the system is governed, together with the metric and the 
EOS, by the Einstein's equation R^'' - ^S.^'^R = nTi'' and 
the equation of motion T'^i-k = [2^ . 
For adiabatic motion, we have. 



d{wui 
ds 



dzu 



.k„,l 



2 dx-" 
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dx^ 



da 
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(4) 
(5) 



where signifies the four velocity, vu and a are the en- 
thalpy and entropy, respectively. For an isentropic pro- 
cess, i.e. a — constant, the term goes to zero. 

The above two equations are the starting point for de- 
riving the appropriate continuity and Euler's equations. 
If we want the equation for we only need to con- 
sider the variation of the gravitational potentials with 



r and the terms containing partial derivatives w.r.t 9 
can be neglected. For our metric, \/—g comes out to 
be y/—g — e^^'^°'r'^ sin9\ the four vectors are also cal- 
culated from the metric. To derive the hydrodynamical 
equations, we use 'auxiliary' time dr instead of ordinary 

^^dt. 



/gtt 



time dt, defined as dr 

Calculating these values from the metric and putting 
them in eqn. (4) and (5), we get the continuity and the 
Euler's equations, given by 



V3 dr 



de I dv 



dv ^ 2v V 

■V — ) -\ h -cot9 

OT r r 

or or 

I ,dp dp. 1 .dv dv. 
ZD dr dr dr dr 



(6) 



where, v is the r.m.s velocity and W is the inverse of 
the Lorentz factor. This system of equations differs from 
its SR counterpart in the appearance of the gravitational 
potentials, describing the gravitational forces acting on 
the fluid in its own gravity field. Once the two basic 
equation are set, we now proceed exactly as in We 
define v as the front velocity in the nuclear matter rest 
frame and n — ^ 

Of- 

speed in the medium, 
relation 

dr „ 

(8) 



■gl is the square of the effective sound 
dr and dr are connected by the 



— ^vG 



dr 



where G is given by 



G 



e7+p _ e''-Pr'^uj^sin'^6 



(9) 



The other parameters of the equation are defined as 
vujrsin9 „ 2uj'^rsin9 2uP'e''~''rsin9 

^ = ^ = ^^ + Ai ' 

Bl vujrsin9 2uje^ Pr'^sin9 vrsin9 

B = —r, — ; C — - 



Al CI 



Al 



CI 



where 



Al = e^+P - e^-Pr^uj^sin^9; 
Bl = -e'^+P - e''-Pr^uj^sin^9] 
CI — \/ r'^ixP'sin?9 — e^P\ 

After a bit of algebra, we get a single differential equa- 
tion for v: 



W^v[K + Kl] 



dv 



dr 2[i;2(l + G)2-n(H-w2G)2] 



(10) 



where, 



r. 2^\/f^7 da 2 cot9 . 

K = 2nil + v^G)i^^+- + - + —) 

Kl^il + G)iA'^+B'/ + c'^ + E) 
dr dr dr 
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Lu = and sin9 = 1 in this equation yield the equation 
for the static star and if we put all potentials equal to 
zero, we recover our equation for the SR case p^ . 

Integrating eq. (10) over r from the center to the sur- 
face, we obtain the propagation velocity of the front along 
the radial direction. The initial condition, i.e, the veloc- 
ity at the center of the star, should be zero from symme- 
try considerations. On the other hand, the 1/r depen- 
dence of the dv/dr in eqn. (10) suggests a steep rise in 
velocity near the center of the star. 

Our calculation proceeds as follows. Having con- 
structed the density profile of the star for a fixed central 
density, eqns. (2) and (3) specify the respective flow ve- 
locities vi and V2 of the nuclear and quark matter in the 
rest frame of the front, at a radius infinitesimally close 
to the center of the star. This would give us the initial 
velocity of the front (— wi), at that radius, in the nuclear 
matter rest frame. With this, we integrate eqn. (10) 
outwards along the radius of the star. The solution gives 
the variation of the velocity with position as a function 
of the time of arrival of the front. Using this velocity 
profile, we can calculate the time required to convert the 
whole star using the relation ^ — vG. 
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FIG. 1: Variation of velocity of the front along the radial 
direction of the star for three different cases, namely SR, non- 
rotating GR and rotating GR. 
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FIG. 2: Variation of velocity of the front along the radial 
direction for different rotational velocity of the star. 



ken curve for non-rotating GR case and the dotted curve 
for the rotating GR case with % = 0, i.e at the equator. 
As expected, for all the three cases, the velocity shoots 
up near the center of the star and saturates at a certain 
velocity for large radii. This is due to the asymptotic 
behaviour of the differential equations. It can also be 
clearly seen that the GR effect increases the velocity of 
the front considerably (maximum by 30%) and the effect 
is most pronounced for the static case. The rotational 
effect of the star seems to suppress the GR effect and 
therefore the velocity of the front decreases. The result 
becomes clearer if we look at fig. (2) where we have plot- 
ted the front velocity with equatorial radius for different 
rotational velocities; as the rotational velocity increases, 
the velocity of the front decreases. 

From fig. (3), we find that the front velocity is max- 
imum along the polar direction and minimum along the 
equator. This is quite an interesting result as the equator 
has the maximum length and the front takes the maxi- 
mum time to convert nuclear matter to two-flavour quark 
matter in that direction than towards the pole. There- 
fore, at any particular instant of time, we may have a 
situation where the polar part of the star has been con- 



For a static star, being spherically symmetric, the 
problem is rather simple; for a rotating star, however, 
the asymmetry has to be taken care of. For that, we 
would have to introduce a new parameter x — cosO, along 
the vertical axis of the star. We start our calculation by 
choosing the central density of the star to be 7 times the 
nuclear matter density, for which the Keplerian velocity 
of the star is 0.67 x 10~'*sec~^ (the rotational velocities 
given in fig. (2) are all in units of lO^^'sec^^). For this 
central density, the initial velocity of the front comes out 
to be 0.45. For a complete understanding of the GR ef- 
fect, comparison with respect to SR is desirable. After 
solving the differential equation for static and rotating 
star, we have plotted in fig. (1) the propagation velocity 
of the front along the radial direction of the star for three 
cases. The unbroken curve is for the SR case, the bro- 
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FIG. 3: Variation of velocity of the front along the radial 
direction for different x- 
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FIG. 4: Variation of time of arrival of the front at certain 
radial distance for different cases. 

verted while along the equatorial direction, the front is 
still propagating. The result is interesting but not unex- 
pected as at the pole the EOS is much steeper than at 
the equator (and we have used the slope of the EOS, 'n' 
to calciilatc the front velocity). At some distance from 
the center of the star, the propagation front breaks up 
into several distinct fronts propagating with different ve- 
locities in different directions. 

From fig. (4) we can see that the time taken by the con- 
version front to convert the neutron star to two-flavour 



quark star is of the order of few ms. The static star 
takes the minimum time (3.3tos) whereas the rotating 
star takes the maximum time (5.1ms) due to the enlarged 
equatorial radius. The polar part of the star needs much 
lesser time for conversion (3.1ms), even less than static 
star, as its radius gets compressed. 

To summarize, we have shown in this letter that the 
conversion of nuclear matter to quark matter in com- 
pact stars, especially rotating stars which arc more real- 
istic than static stars, is strongly affected by GR effects. 
The emergence of different conversion fronts, propagating 
with different velocities along different radial directions, 
is a novel finding, not anticipated by Newtonian or SR 
calculations. It remains to be explored whether the in- 
corporation of dissipative effects materially changes the 
conclusion. Intuition would point to a dampening of mo- 
tion under dissipation, but the relative gradients between 
the different fronts may also be accentuated, with impor- 
tant physical consequences. Involved though it is, such 
an investigation is on our immediate agenda. 
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